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LONG-TERM  GOALS 

Out  interest  is  in  understanding  how  stratification,  turbulence,  and  circulation  interact  in  estuarine 
flows.  It  is  our  hypothesis  that  the  role  of  stratification  in  suppressing  turbulence  is  central  to  this 
interaction.  Thus,  a  major  focus  of  our  work  is  in  finding  simple  relationships  (if  possible)  connecting 
local  flow  conditions  to  turbulence  state  and  turbulent  fluxes.  We  are  also  interested  in  using  field  data 
from  complex  estuarine  flows  to  test  circulation  models,  both  in  terms  of  straightforward 
comparisons/calibration  and  in  terms  of  improving  the  ability  of  circulation  models  to  properly 
represent  the  physics  of  these  flows. 

OBJECTIVES 

This  project  is  focused  on  using  data  collected  during  the  July  2006  COHSTREX  experiment  in  the 
Snohomish  River  Estuary  (SRE)  to  address  the  following  issues: 

(1)  Construct  the  momentum  balance  through  the  tidal  cycle  and  including  spring  neap  variations  in 
tidal  forcing  as  well  as  an  assessment  of  the  uncertainty  of  each  term. 

(2)  Analyze  the  stratified  turbulence  data  acquired  with  the  ADCPs  in  the  SRE  to  determine  when 
and  where  mixing  of  salt  and  momentum  occur  in  the  system  and  thius  to  understand  the  role  of 
the  interaction  of  stratification  and  turbulence  in  shaping  flows  in  the  SRE. 

(3)  Construct  an  analytical  framework  for  studying  residual  circulation  in  macro-tidal  estuaries, 
focusing  on  simplified  geometries  and  forcing  (e.g.  constant  baroclinic  pressure  gradient),  and 
use  that  framework  to  examine  our  field  data. 

(4)  Understand  the  dynamics  of  what  we  suppose  to  be  trapping-driven  frontogenesis  in  the  complex 
near-field  region  of  the  SRE  that  was  the  focus  of  the  2006  COHSTREX  experiment. 

APPROACH 

During  EY2010  our  efforts  have  been  focused  on  further  analysis  of  the  data  obtained  in  the  July  2006 
COHSTREX  experiment  in  the  Snohomish  River  Estuary  in  Washington.  This  entails  using  the  data  to 
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construct  budgets  for  salt  and  momentum  that  will  enable  us  to  understand  where  and  when  particular 
physical  processes  like  tidal  straining  are  important.  We  also  have  been  foeusing  on  analysis  of  our 
data  in  terms  of  tidal  and  intra-tidal  variability  of  the  flow  when  viewed  in  sigma  eoordinates,  i.e.,  in 
eoordinates  that  are  defined  as  a  fraetion  of  the  instantaneous  depth,  an  analytieal  framework  made 
neeessary  by  the  large  variations  in  depth  observed  in  the  Snohomish.  This  work  forms  the  thesis  of 
Ms.  Sarah  Giddings  whose  thesis  is  eo-supervised  by  the  PI  and  Dr.  Derek  Fong. 

WORK  COMPLETED 

Our  aetive  areas  of  researeh  in  FY2010  have  ineluded  detailed  analysis  of  both  the  intratidal  variability 
and  the  subtidal  dynamies  in  the  Snohomish  River  Estuary  (SRE).  On  the  tidal  time  seale  we  foeused 
on  the  evolution  of  stratifieation,  vertieal  mixing,  and  longitudinal  dispersion,  as  well  as  frontal 
dynamics.  This  analysis  provided  good  insight  to  the  dynamies  in  similar  strongly  foreed,  strongly 
stratified  estuarine  systems  whieh  are  prevalent  yet  less  well  studied  than  partially  mixed  systems.  We 
also  have  been  working  on  investigating  the  subtidal  variability  of  the  SRE,  aiming  to  extend  our 
understanding  of  estuarine  dynamics  to  the  less  well  understood  strongly  foreed,  strongly  stratified, 
large  amplitude  regime  where  eommon  approximations  sueh  as  a  eonstant  horizontal  salinity  gradient, 
dS/dx,  or  a  small  amplitude,  77//i«l  do  not  hold. 

The  work  eompleted  in  the  past  year  has  been  ineorporated  into  one  submitted  manuscript, 
presentations  at  several  eonferenees,  and  two  papers  in  preparation. 

RESULTS 

Intratidal  dynamics  -  stratification  and  vertical  mixing 

Due  to  the  salt-wedge  nature  of  the  system,  both  along-stream  adveetion  and  eoneomitant  straining  of 
the  density  field  dominate  temporal  and  spatial  variations  in  stratification  (Eigure  1).  The  dynamies  of 
the  straining  and  adveetion  vary  substantially  on  both  a  tidal  and  a  fortnightly  time  seale  leading  to 
strong  temporal  variability  in  stratifieation,  vertieal  mixing,  and  longitudinal  dispersion  that  differs  in 
important  ways  from  observations  in  partially-mixed  estuaries.  In  partieular,  flood  tide  adveetion  is 
counteracted  by  straining  and  vertieal  mixing  sueh  that  stratifieation,  shear,  and  longitudinal  dispersion 
are  minimal.  In  eontrast,  ebb-tide  straining  and  adveetion  enhanee  vertieal  stratifieation  and  vertieal 
shear  whieh  enhanees  longitudinal  dispersion  and  leads  to  strong  ebb-tide  interfaeial  mixing. 

Intratidal  dynamics  -  Frontogenesis 

Complex  bathymetry  in  the  SRE  leads  to  trapping  of  mid-density  water  over  intertidal  mudflats  whieh 
later  eonverges  with  dense  water  in  the  main  channel  forming  a  sharp  front.  The  frontal  density 
interfaee  is  maintained  via  eonvergent  transverse  eirculations  driven  by  both  lateral  baroelinic  and 
eentrifugal  foreing.  Importantly,  the  front  leads  to  enhaneed  stratifieation  and  suppressed  vertieal 
mixing  at  the  end  of  the  large  flood  tide.  Despite  these  loeal  effeets,  this  front  does  not  signifieantly 
alter  longitudinal  dispersion.  This  frontogenesis  meehanism  fits  within  the  broader  context  of 
bathymetrieally  driven  meehanisms  in  whieh  varying  bathymetry  drives  lateral  eonvergenee  and 
baroelinie  foreing  and  we  expeet  it  to  oeeur  in  other  systems,  potentially  even  eontributing  to 
longitudinal  dispersion  in  estuaries  with  larger  trapping  regions.  In  this  investigation,  we  employed  the 
entire  suite  of  COHSTREX  measurements  including  numerical  simulations,  RiverRad,  TIR  imagery. 
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REMUS  and  our  other  in-situ  instruments.  The  manuseript  ineorporating  the  results  of  this 
investigation  involved  eollaboration  with  several  other  COHSTREX  investigators  whom  are  ineluded 
as  authors. 


Figure  1:  Stratification  Evolution  at  mooring  MSB:  mean  velocity  and  water  level  (a,  e),  top  bottom 
stratification  (b,f),  and  the  depth-average  of  terms  contributing  to  the  evolution  of  stratification 
(c,  g  and  d,  h)  for  a  spring  tide  (left  panels)  and  neap  tide  (right panels). 

Subtidal  dynamics 

The  most  reeent  area  of  researeh  has  been  the  investigation  of  subtidal  dynamies.  Most  estuarine 
systems  exhibit  residual  (subtidal)  eirculation  driven  by  the  baroelinic  foreing,  i.e.  upstream  at  depth 
and  downstream  at  the  surfaee  whieh  is  eonsidered  a  eritieal  faetor  in  residual  transport  proeesses. 
Residual  eireulation  is  diffieult  to  interpret  however  in  a  shallow,  maerotidal  system  where  the  tidal 
amplitude  is  comparable  to  the  mean  depth.  Profiles  of  residual  circulation  and  residual  density 
structure  are  shown  in  Eigure  2  below.  We  have  been  working  on  several  approaches  including  data 
manipulation  and  analytical  work  to  try  to  interpret  the  observed  residual  circulation  structure, 
compare  it  to  traditional  theories,  and  the  intratidal  processes  discussed  above  that  contribute  to  the 
residual. 

The  residual  circulation  does  not  exhibit  the  traditional  estuarine  structure  (Figure  2a),  rather  there  is 
out-estuary  flow  at  the  bed  and  peaks  mid  water  column.  If  we  convert  our  coordinate  system  to  a 
normalized  depth  coordinate  system,  o=zlD  where  D  is  the  temporally  varying  total  water  depth,  the 
residual  circulation  profile  cleans  up  substantially  (Figure  2c)  and  in  fact  resembles  more  of  a 
traditional  estuarine  circulation.  Similarly,  the  density  profiles  exhibit  the  more  expected  residual 
profiles  (Figure  2d)  as  opposed  to  the  z  coordinate  system  where  they  actually  exhibit  unstable 
residual  profiles  (Figure  2b). 
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Figure  2:  Residual  (tidally  averaged)  circulation  (left  panels)  and  density  profiles  (right  panels)  in 
traditional  z  coordinates  (upper  panels)  and  depth  normalized,  [J coordinates  (lower  panels). 

The  normalized  depth  eoordinate  system  is  useful  for  this  type  of  estuary  beeause  it  eliminates  large 
amounts  of  missing  data  when  the  water  level  is  beneath  high  water.  We  performed  an  empirieal 
orthogonal  function  (EOF)  analysis  on  the  data  in  cr  coordinates  and  found  that  the  first  two  modes 
represented  over  99%  of  the  variability  in  the  data  and  correlated  significantly  with  the  depth  averaged 
velocity  and  the  vertical  stratification  respectively  suggesting  that  EOFl  represents  a  baro tropic  mode 
while  EOF2  represents  a  baroclinic  mode  (Figure  3).  We  note  that  EOF2  also  includes  some  shear 
driven  by  mechanisms  other  than  baroclinic  forcing,  therefore  we  will  refer  to  it  as  a  shear  mode. 

We  also  included  the  Stokes  wave  transport  velocity  with  the  Eulerian  residuals  (4)  to  approximate  the 
Fagrangian  residual.  We  find  that  the  Stokes  wave  transport  is  directed  upriver  at  all  depths  such  that 
the  Fagrangian  two-layer  circulation  (not  shown)  is  stronger  than  it  appears  in  the  Eulerian  residuals. 

To  try  to  compare  the  observed  residuals  to  traditional  estuarine  theory,  we  have  examined  analytical 
solutions  to  the  tidally  averaged  momentum  equations  in  cr  coordinates.  Although  the  crprofiles 
resemble  traditional  theory  we  note  that  the  approximations  required  to  derive  the  results  are  often  not 
met  in  the  SRE,  particularly  near  the  water  surface.  Additionally,  we  have  examined  the  tidally  varying 
momentum  equations  in  detail  and  hypothesize  that  like  many  estuarine  systems,  this  residual  is  driven 
by  intermittent  pulses  of  two  layer  circulation  created  by  straining,  intratidal  variability  in  vertical 
mixing,  and  the  strongly  temporally  varying  horizontal  salinity  gradient  that  we  investigated  in  our 
submitted  manuscript.  Overall  it  appears  that  this  cr  coordinate  approach  shows  promise  for 
interpreting  residual  circulation  in  systems  with  a  large  tidal  amplitude  and  linking  these  systems  to 
traditional  estuarine  theory. 
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Figure  3:  The  upper  panels  show  the  depth  averaged  along-stream  velocity  (black  line)  along  with 
the  amplitude  of  EOFl  (dotted  black  line).  The  lower  panels  show  the  vertical  stratification  (black 
line),  horizontal  density  gradient  multiplied  (grey  line)  along  with  the  amplitude  of  E0F2  (dashed 
black  line).  The  left  panels  show  spring  tides  while  the  right  shows  neaps. 
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Figure  4.  EOF  components  of  the  residual  circulation 


IMPACT/APPLICATIONS 

The  Navy  is  responsible  for  conducting  operations  in  complex  coastal  environments,  such  as 
macro  tidal  systems  found  on  the  west  coast  of  the  USA  or  on  the  coasts  of  eastern  Africa,  northwest 
Australia  or  Korea.  The  work  in  this  project  will  help  build  the  dynamical  framework  for  these  flows 
that  is  needed  to  model  them  or  to  interpret  data  acquired  (e.g.  remotely)  about  them. 
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This  project  is  a  continuation  of  the  ONR  MURI  COHSTREX  (N00014-05-1-0485)  involving 
researchers  from  the  University  of  Washington  and  Stanford. 
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